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The renormalization group is applied to the <j> 4 model in the symmetry broken 
phase in order to identify different scaling regimes. The new scaling laws reflect 
nonuniversal behavior at the phase transition. The extension of the analysis to 
finite temperature is briefly outlined. It is mentioned that the coupling constants 
can be found in the mixed phase by taking into account the saddle points of the 
blocking procedure. 



1 Introduction 

In the search of the quark confinement mechanism the haaron model has been 
proposed! because it comprises the lesson to be learned from lattice QCD. The 
characteristic feature of the model is that the linearly rising potential between 
static color charges arises from a simple sine-Gordon type effective model after 
a partial resummation of the Haar-measure vertices of the path integral. This 
was a rather puzzling result since it was difficult to accept that the leading long 
range force comes from vertices which are nonrenormalizable and irrelevant. 
But one can easily find the explication of this apparent paradox: On the one 
hand, the renormalizability stands for the relevant or marginal behavior of the 
operators in the ultraviolet scaling regime. On the other hand, the confining 
forces are observed beyond this scaling regime, where new scaling laws arise at 
the low energies. There is no reason to expect that the relevant operator set 
agree for the ultraviolet and the infrared scaling regime. Thus one is left with 
a more general question whether the set of relevant operators in a model may 
differ in different scaling regimes and what consequences such a phenomenon 
might have. 

When the mass scale is explicitly given in the lagrangian as in the massive 
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4> A model then the simple perturbation expansion is sufficient to study the 
different scaling regimes. There are cases when a partial resummation of the 
perturbation expansion is needed to generate the mass gap, like for photons 
at finite temperature. In these cases where perturbation expansion applies 
there is no transmutation of the degrees of freedom, i.e. one finds the same 
particles at every energy and each fixed point is Gaussian. When the low energy 
scaling law is accompanied by the appearance of new particles, such as in the 
two dimensional Gross-Neveu model or in QCD then new relevant operators 
are expected which are responsible for the formation of the new composite 
particles and the low energy fixed point is not Gaussian anymore. The generic 
mechanism for the nonperturbative modification of the scaling law and the 
generation of new relevant operators is the condensation. This phenomenon 
may occur either at the Iowa or at the high energjH scaling regimes. Since 
the mass generation is usually achieved in High Energy Physics by the help 
of the condensates one may find similar complication in the unified models, 
as well, despite their perturbative appearance. The present contribution is a 
brief summary of some results obtained in this direction. 

We shall first argue in Section 2 that the appearance of several fixed points 
and scaling regimes is a rule rather than an exception in High Energy Physics. 
The powerful Wegner-Haughton form of the renormalization group equation is 
introduced in Section 3 as the method to tackle our problem in the symmetry 
broken phase of the </> 4 model. A diverging and a focusing effect of the renor- 
malization group flow is discussed in Section 4. Section 5 is a brief digression 
into the structure of the mixed phase. The generalization of our results to finite 
temperature is the subject of Section 6. Finally Section 7 is for the summary. 

2 Multiple Fixed Points 

The models with intrinsic mass scales possess at least two distinct scaling 
regions, one at the ultraviolet and another one at the infrared side of the mass 
scale. The infrared scaling is usually called trivial because it can be proven 
that there is no non-Gaussian relevant or marginal operator. In fact, for models 
with finite correlation length the fluctuations are exponentially suppressed at 
large distances and the evolution of the running coupling constants slows down 
in the infrared limit. The manifold of the possible attractive infrared fixed 
points is parametrized by the initial values of the relevant coupling constants 
of the ultraviolet scaling regime given at the scale of the ultraviolet cutoff. 
When the excitation spectrum has no gap above the vacuum or there is an 
instability then the long range interactions might be so strong as to change 
this simple situation. The result is that divergences might pile up and drive 
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Figure 1: The renormalized trajectory of the Theory of Everything (TOE) starts at the 
(supposed) ultraviolet fixed point at k = A, it passes by the fixed points of the Grand Unified 
Models (GUT), the unified Electro- Weak theory (EW), the strong interactions (QCD), the 
electromagnetic interactions (QED) certain fixed points of the Solid State and Condensed 
Matter Physics (CM) and finally approaches the ultimate IR fixed point, k — > 0. The circles 
denote the domains of lincarizability, the asymptotic scaling regions. 

the trajectory away to another, unexpected region in the coupling constant 
space. These runaway trajectories indicate the existence of a relevant operator 
in the infrared scaling regime. 

The presence of several scaling regimes is easily recognizable at the The- 
ory of Everything. Whatever theory will proven to be that, its renormalized 
trajectory should be traced down in a space which contains all coupling con- 
stants what is used in physics at finite energies. From the coupling constants 
of possible composite models at so far unexplored high energies through the 
parameters of the Standard Model down to the coupling constants in Solid 
State Physics one includes several axes in this space. On the renormalized 
trajectory depicted in fig. 1 one observes the scaling laws characteristic of 
different interactions in the energy regime where the trajectory is in the lin- 
earizability region of a fixed point. Note that the trajectory may be influenced 
by the environment and bifurcates into different thermodynamical phases in 
the IR regime in different laboratories. 

It is mathematically certainly correct to say that the renormalized tra- 
jectory, the set of the physical "constants", is given by the initial condition 
at the TOE. In this respect the High Energy physicist seeks the few ultimate 
constants of Nature. But the same algebra of observables is classified at each 
scaling region according to the actual scaling laws and there is no guarantee 
that the set of relevant and marginal operators is found to be always the same. 
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The renormalizable (relevant or marginal) operators of a scaling regime are 
usually nonrenormalizable (irrelevant) at the higher energy fixed points such 
as the quark-gluon QCD vertex appears as a nonrenormalizable one in a com- 
posite model for the quarks and gluons. Consider now a coupling constant, 
denoted by g n - r (k) what is nonrenormalizable at high energy and becomes 
relevant at a lower energy scaling regime. Then this coupling constant under- 
goes a suppression at the high energy scaling regime. How can it became a 
relevant coupling constant at the low energy scaling regime? What is its role 
during the evolution? According to the usual scenario the initial value of the 
irrelevant coupling constant, g n - r (A), modifies the theory in an overall scale 
and the relation between observables of the same dimension is given by the 
initial value of the renormalizable coupling constants only. This is expressed 
by the condition for the beta function for any coupling constant g as 

hm - — 2—— = lim k—— '-— = 0, (1) 

A/fe^oo dg n - r (A) A/fc^oo okog n - r (A) 

where the coupling constants are made dimensionless by the help of the running 
cutoff, k. 

One can imagine the following two, opposite possibilities as far as the 
coupling constant g n -r is concerned: 

1. Divergence: The universality, eq. ([!]) is violated because the trajecto- 
ries with slightly differing initial conditions for g n - r diverge from each 
other. The initial value of g n - r must then be specified at the TOE and 
it becomes an independent free parameter. 

2. Focusing: The dimensionless quantities at low energy are (at least locally) 
independent of the initial values of the coupling constants of the TOE 
and are determined by one of the lower energy fixed points. This is 
a strong version of the universality because the coupling constants at 
low energy are independent of the initial values of the renormalizable 
coupling constants. 

The existence of different scaling regimes may lead to serious problems in 
indentifying the important parameters of the theory In both cases mentioned 
above the determination of the physical content of the theory at different en- 
ergies in terms of the initial values of the relevant coupling constants of the 
TOE is, though being mathematically possible, unfeasible by means of mea- 
surements with small but finite errors. 

It is instructive to consider a simpler model with two scaling regimes whose 
generic scaling patterns are listed in table 1. Consider the theory of photons, 
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Table 1: The four classes of the coupling constants in QED. 
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Figure 2: The qualitative dependence of the running coupling constants in the function of 
the cutoff, a = 2tt / A. The UV and the IR scaling regimes are shown. The coupling constant 
is supposed to be constant in between for simplicity. See Table 1 for the details. 
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electrons and nuclei in the presence of a chemical potential for the baryon num- 
ber. For certain values of the chemical potential the ground state is a solid 
state lattice. We identify two scaling regimes, the ultraviolet one characteris- 
tic of QED in the trivial, homogeneous vacuum and an infrared one which is 
governed by the lattice effects. The electron mass, m e , is a relevant parameter 
in each scaling regimes since it is a renormalizable coupling constant of QED 
and appears in the equations of Solid State Physics. The muon mass, m M , is as 
renormalizable parameter at high energy as m e but drops out from the physics 
of the solids because the processes (without neutrinos) at the energies of eV 
the muon contributions are always dominated by the contributions from the 
electron. The six-electron vertex is irrelevant everywhere. The four electron 
vertex with the coupling constant G4 is the most interesting interaction. It is 
a nonrenormalizable, irrelevant vertex of QED. But it becomes relevant at low 
energies since it is the effective vertex which is generated by the attractive force 
between the Cooper-pairs and drives the transition into the superconducting 
state. In the simple perturbative treatment of QED G4 is suppressed in the 
ultraviolet scaling regime so its high energy initial value in the QED lagrangian 
is set to zero. It was only after the experimental discovery of the superconduct- 
ing phase and its explication by the BCS ground state that the importance of 
G4 at lower energies was demonstrated by means of the partial resummation 
of the perturbation expansion. This is an example where the renormalization 
group was used to find a coupling constant what appears to be unnecessary 
at high energies but nevertheless is important at low energy. Such parameters 
will be called hidden coupling constants after their undetectable small values 
at intermediate energies. 

The divergence, case 1 mentioned above, requires that the suppression of 
g n - r in the ultraviolet regime is weaker than the amplification during the low 
energy scaling. If the fixed point of the infrared scaling has a certain region 
of the ultraviolet coupling constants in its attractive zone then focusing, the 
case 2, is realized. We present here a study of the four dimensional </> model 
in search for the manifestations of these phenomena. 

3 Wegner-Haughton Equations 

In order to study the role of nonrenormalizable operators on the evolution one 
has to be able to handle the mixing of a large number of operators during the 
renormalization group transformation. This is achieved in an elegant manner 
by the Wegner-Haughton equation This is a functional differential equa- 
tion describing the evolution of the renormalized, blocked bare action under 
the change of the cutoff. According to the usual strategy the running cou- 
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pling constants are identified with the bare ones and the observational scale 
with the cutoff so the evolution of the bare coupling constants qualitatively 
reproduces the trajectory for the running coupling constants. We shall be 
satisfied here to indicate the derivation of the leading order equation in the 
gradient expansions, the preliminary results indicate that the higher orders do 
not change our conclusions. 

Let us write the action corresponding to the ultraviolet cutoff k as 



S k [<f>] = J d d x 



(2) 



2 

According to the usual Wilson-Kadanoff blocking procedural we write 

to'M = / D[4>']e-^+^ (3) 



e 



where k! < k in the Euclidean space-time. The Fourier transform of the field 
variable <f> and cf>' are nonvanishing for p < k' and k' < p < k, respectively. 
The right hand side is evaluated by means of the loop expansion, 

Svffl = S k [cj ) + ^ } + ^trlog5 2 S + O(h 2 ), (4) 

where 

and the saddle point, <f>' , is defined by the extremum condition 

H'{p) ' [ ' 

in which the infrared background field, 4>(x), is held fixed. Eq. (0) is the 
generalization of the Wegner-Haughton equation for condensates. One can 
prove that the saddle point is trivial, </>(, = 0, so long as the matrix <5 2 5(x, y) is 
invertible and the infrared background field is homogeneous, <fi(x) = $. What 
is remarkable in this equation is that each successive loop integral brings a 
suppression factor 

k d ~k' d „{k-k'\ 

due to the integration volume in the momentum space. By choosing an in- 
finitesimal fraction of the degrees of freedom to be eliminated in a step we 
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find a new small parameter, 6k/ k' = k — fc'/fc', suppressing the higher loop 
contributions in the blocking relation. 

In the leading order of the gradient expansion, the so called local potential 
approximation, the only function characterizing the action is the local poten- 
tial, Uk{4>) : so we choose the infrared background field homogeneous, </>(x) = 
and obtain from (||) 



£/ fc _ 5fc ($) = E/ fe ($) + - 



where we introduced the notation 



d d p 



(27T) 



d 2 U k (<i>) 

■ 



In the limit 5k — > one ends up with the differential equation 
4^) = -^lo S [^ + W)] 



dk 



2(2tt) 



with 



2d JT n/2 ' 



(8) 



(9) 



(10) 



(11) 



the projectiorll of the functional equation (Q) onto the functional form (||). It 
is instructive to expand this equation in — ra 2 ., where 



CT"(0), 



when we find 



log(fc 2 +m£) + 



i ( u'jm-™\ 

k 2 + m? 



(12) 



(13) 



One recovers here the usual one loop resummation of the effective potentials 
except that the loop momentum is now restricted into the subspace of the 
modes to be eliminated. Note that the derivation of eq. ( |l0| ) shows that the 
restoring force for the fluctuations into the equilibrium is proportional to the 
argument of the logarithm function. Thus when 



fc 2 + U% ($) = 

then nontrivial saddle point should be used. The original equation, 
always valid because the action is bounded from below. 



(14) 
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We define the coupling constant and the beta function to the n-th order 
vertex for the vacuum 4>{x) = $ as 

ft = k ^ n (k) = ^^ k u k m, 

where in the last equation we assumed the analycity of the potential in k and $, 
what holds except at the singular points. By taking the successive derivatives 
of the renormalization group equation (^) one obtains 



--^ 1 V n (G 2 ,---,G n+2 ), (Hi) 



where 

g_ = 

k 2 + g 2 



Gn = T^hr (17) 



and 

is a polynom of order n/2 in the variables Gj, j = 2, ■ ■ ■ ,n + 2, 



ft = 








ft = 




r 2 




ft = 


G 5 


— 3G3G4 


+ 2GI 


ft = 


^6 


- 4G 5 G 3 




ft - 


G 7 


— 5G6G3 


- 4G 5 G 4 + 20G§G B - 6G^G 5 + 30G^G 3 



(19) 



-60G^G 4 + 24G|, 



etc. It is interesting to verify that V n contains the integrand of all one loop 
graphs what contribute to (3 n . 

We change now to dimensionless parameters, 

k->j, U^k d U, <j>^ki- l <t>, gn^k^-i^gn, (20) 

and 

G„ — T ^-, (21) 



1 + 92 



what will be used below. 



One can distinguish an ultraviolet and an infrared scaling regime, for 
k 2 » and for k 2 « |m||, respectively. In the former one recovers the 
usual renormalization group coefficient functions used in studying the asymp- 
totic scaling. The latter is trivial for Wq > when the factor k d suppresses the 
evolution at the infrared fixed point. The infrared scaling law is presumable 
trpial in the case m 2 (0) = because the vacuum of this model is supposed to 
bco at $ 7^ where the fluctuations become massive again. 

The symmetry broken phase^ is characterized by the condition that there 
is a scale, fc cr (<!>), where the spinodial instability occurs and the restoring force 
of the fluctuations is vanishing, 

k 2 cr (<S>) = ~U'>J<P) (22) 

in a region around $ = 0. We shall call the line k 2 r {&) on the plane ($,/c 2 ) 
critical. The renormalized trajectory has discontinuous derivatives along the 
curvcl and the system is in a mixed phase for k 2 < k 2 r (Q). The tree level 
instability induces nontrivial saddle points for the blocking ^ and eq. (|To|), 
being based on the vanishing of the saddle point, is no longer valid. As the 
critical curve is approached in decreasing k during the blocking the denomina- 
tor 1 + g2 (k) becomes small in the beta functions and new scaling laws develop 
as a precursor of the mixed phase. There is no reason to expect that these new 
scaling laws should be trivial. We shall study the renormalization group flow 
in the vicinity and below the critical line. 

One can nowadays experience a renaissance of the infinitesimaL renormal- 
ization group methods. The so called exact renormalization groupLH is similar 
in spirit to the Wegner-Haughton equation but it follows the evolution of the 
generator functional for the connected or the 1PI vertices. The advantage of 
this method is that it produces directly the particle physics motivated run- 
ning coupling constants which are based on the scattering amplitudes. The 
bare renormalization group, (j^), yields simpler expressions for the evolution of 
the bare coupling constants of the action. We found this latter method more 
attractive since, as mentioned at the beginning of this section, the evolution 
of the bare and the running, renormalized coupling constants is qualitatively 
similar. 

Another distinguishing feature of the blocking procedure employed in this 
work is that the momentum space cutoff is imposed in a sharp manner. The 
smooth cutoff is believed to be superior and its use is more widespread. We 
shall argue below that the problems with the sharp cutoff are not incurable with 
careful methods and that actually no smooth cutoff is applicable for models 

a It could naively be identified with m 2 (0) < but we should not forget the Maxwell 
construction what sets m 2 (0) = 0. 
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with instabilities. The sharp cutoff induces diffraction integrals during the 
blocking what represent oscillating forces at large distances. The oscillations 
cast doubt on the physical significance of the running parameters of the blocked 
action and were blamed for the occurrence of the infrared singularities in the 
renormalized trajectory. Attempts to eliminate the oscillating part from the 
blocked action lead to the introduction of a smooth cutoff. Let us deal with 
the cases where the singularity occurs at k 2 cr = and k 2 cr > separately and 
consider a physical quantity obtained in the framework of the loop expansion, 



where I n {e, k) stands for the n-th order loop integral with the range of inte- 
gration e < \pj\ < k, j = 1, • • • , n. In the case of the infrared unstable theories 
it is necessary to introduce the infrared cutoff, e, what is removed after the 
computation is completed. The blocking transformation with sharp cutoff can 
be used to obtain the right hand side of the differential equation 



The integration of this equation yields the quantity sought when e — > 0. So 
long as the nontrivial vacuum, e.g. the condensate, what shields the infrared 
divergences is properly incorporated in the computation the thermodynamical 
limit is well defined and the e-dependence is continuous at e = 0. Moreover 
any infrared cutoff should yield the same thermodynamical limit. When the 
singularity occurs at finite scale then the quantity P(0, k) can be thought as it 
had been obtained in the effective theory with a sharp ultraviolet cutoff at A = 
k. The bare parameters of this effective theory have singular cutoff dependence 
at k = k cr . One might argue that this singularity is an artifact of the cutoff 
employed because the observables in the effective theory, being renormalization 
group invariant, show no singular fc-dcpcndcncc. But we see no conceptual 
problem with singular, i.e. nondiffcrcntiable renormalized trajectories if the 
singularity corresponds to a real physical effect, i.e. some instability what 
shows up at a well defined, sharp value of the momentum. One should simply 
make sure that the physical effects behind the singularity have properly been 
accounted for during the solution of the effective theory or the continuation of 
the blocking procedure. 

The argument prohibiting the application of any smooth cutoff in theories 
with condensate or other instabilities is the following. In order to show the 
exactness of the infinitesimal renormalization schemes one has to assume first 



oo 




(23) 



n=l 
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the validity of the loop expansion. In fact, the counting of the power of 5k is 
done in the loop expansion where one has to integrate around the saddle point. 
Without the proper choice of the saddle point the formal steps in arguing about 
the exactness are not valid. The point is that the subsequent elimination of 
the degrees of freedom modifies the integrand for the unstable mode. Due to 
the factor Ti~ l coming from the nontrivial saddle point the loop corrections 
of the stable modes what are computed after the elimination of the unstable 
mode yields the contribution O(h ). Thus all stable modes should completely 
be eliminated before one arrives at the unstable sector of the theory. The 
models with condensate require a blocking method where the stable modes 
are eliminated completely before arriving at the instability and the tree level 
structure of the saddle point expansion must be retained. Finally we note that 
this type of cutoff poses no problem in going to higher orders in the gradient 
expansion. 

4 Zooming into a fixed point 

Our goal is to follow the evolution of the blocked action from the initial condi- 
tion set at k = A — 1 towards the infrared regime. We shall consider the one 
component four dimensional (f> 4 model in the symmetry broken phase. Thus 
the integration of the differential equation (|l^) in k runs into a singularity 
for $ = at k = k cr (0) > 0. When the evolution in the outer, stable re- 
gion is approximated by the tree level expression, £4($) = U\(&), then the 
singular line, fc^ r (<I>), is an upside down parabola, k 2 cr = —m\ — g4(A)$ 2 /2, 
on the plane ($, fc 2 ). But this is an inconsistent approximation for the tree 
level saddle point structure is nontrivial in the unstable region for k < fc cr (<i>). 
Since the input for the elimination of a mode at k is the potential Uk+Sk{^) 
for — co < $ < oo the result of the naive blocking what does not take into 
account the nontrivial saddle point structure is built on the wrong potential 
and is incomplete by a term 0(U) even in the outer, stable region. We shall 
avoid the problem of the singularity by assuming that the potential is 
analytical inside and outside of the unstable region and its only nonanalytical 
behavior is confined on the curve fc 2 r (<I>). This assumption allows us to study 
the evolution equation locally in <i>, i.e. by expanding the potential as 



and following the evolution of the coupling constants g„(<I>o) what obeys the 
loop expansion with the trivial saddle point in the outer, stable region. 




(25) 
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We note that the Taylor expansion of the potential motivated here by 
the decoupling of the stable and the unstable region proved to be necessary 
in the numerical integration of eq. (|l0|). The numerical integration of this 
equation by simply discretizing the variables $ and k becomes highly unstable 
at k 2 « k 2 r (0) because the logarithm function amplifies the numerical errors 
in computing ($). The smoothing or interpolating techniques we tried were 
unable to stabilize the solution. Thus we finally integrated numerically the 
coupled differential equations imposed at $o = 0, 

da k 4 

= -( n ~ 4 )f« " V n(G 2 , ■ ■ ■ , G n+2 ), (26) 

for n = 1,---,N. 

We can present here only some preliminary numerical results, the detailed 
account will be given elsewherai3. They suggest the existence of two distinct 
scaling regimes, an ultraviolet, k 2 >> k 2 r (Q), and another one in the vicinity 
of the singular line. Though the singular line is not a fixed point since k 7^ 
or cxd nevertheless the singularities suggest to parametrize the renormalized 
trajectory by 




This parametrization possess scaling region at at k = and 00. Furthermore we 
shall argue below that the potential is renormalization group invariant in the 
interior, unstable region. In this manner the whole unstable region represents 
a single fixed point. Note that though the critical point belongs to finite scale, 
k = fccr(O) 7^ 0, the singularities presented below require a dense enough 
spectrum for the momentum operator, the execution of the thermodynamical 
limit. 

The result of the numerical integration indicates an attractive fixed point 
at k — fc CI .(0) for N — 10, i.e. by truncating the potential at 0($ 20 ). The value 
of 8k was adjusted during the integration in the range 10~ 18 < Sk < 10~ 2 . 
We compared the third and the fourth order Runge-Kutta approximation for 
the coupling constants and 5k was chosen to keep the relative local error on 
them less than 10~ 15 . After the system leaves the ultraviolet scaling regime 
the higher order coupling constants undergo oscillations with large amplitude 
and a new scaling law is found in the vicinity of the instability. For example 
<?2o reaches the range of 10 -12 for k « k cr (0) after having gone through the 
peaks at 10 23 . The fixed point at k = corresponds to the potential 

U kcr(a) (<Z>) = ~k 2 cr m 2 - (28) 
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The infrared fixed point is trivial and its attractive zone seemed to extend 
over the whole symmetry broken phase of the renormalizable <j> model. The 
approach to the fixed point is such that G n — > 0, n > 2 as k 2 — > k 2 r (0). This is 
consistent with the observation that the fourth equation of (|l9|) excludes any 
finite, nonzero value for (74 at the critical point what belongs to the trivial root 
of V n , n> 2. Thus the scaling at this critical point shows strong universality, 
its result does not depend even on the renormalizable coupling constants. The 
symmetry broken </> 4 model when its potential is truncated at 0(cf> 20 ) realizes 
an example of the case 2 mentioned in Section g. 

When the potential is truncated at 0($ 22 ) or at higher order the qualita- 
tive behavior of the solution is different. At a certain point the trajectory with 
N > 11 suddenly turns away from the solution with N < 10 and the coupling 
constants start to diverge violently towards —00. The qualitative behavior 
seen on fig. 3 remains the same except the singularity becomes stronger for 
N > 11. The sudden departure of the trajectories with N = 10 and 11 can be 
traced back to the contribution of 322 to P20 which happens to be exceedingly 
large at k « k cr (0) and starts to push the lower order coupling constants one 
after the other towards —00. The integration of such a singular curve requires 
extreme numerical accuracy. The quadruple precision was used in the codes 
and the relative precision of the beta functions was kept below 10~ n . The 
values of 5k were between 10 -17 and 10 -13 at k ss fc cr (0). 

It is interesting to observe the derivative of few beta functions with respect 
to the initial value of a nonrenormalizablc coupling constant, 36 (A), 

what is presented in fig. 3. The derivatives are small in the ultraviolet scaling 
regime and diverge at the instability. This divergence shows the presence of 
a strong amplification mechanism due to the instability what lets the small 
details of the microscopic interaction felt in the vicinity of k = k cr (0). By 
assuming the singular behavior kr Vn for (|2^) as far as the k dependence is 
concerned at k « fc cr (0) and the perturbative 0(A 2 ) suppression mechanism 
in the ultraviolet scaling regime for g§, one arrives at 

d 9e (A)~ k A ' (30) 

what reveals a sensitivity of the dynamics at the scale k to the values of the 
nonrenormalizable parameters at 

A«r-' 2 . (31) 
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Figure 3: The evolution of log \g n \ for some values of n at N = 11, with the initial conditions 
32(A) = -0.1A 2 , 94 (A) = 0.01, ge(A) = 10~ 4 A~ 2 and 9j -(A) = 0.0 for j = 4,---, 11. Where 
g n (k) changes sign during the oscillations for n > 4 there log|g n (fc)| shows a cusp whose 
altitude is finite due to the finite resolution of the k values on the plot. The points shown 
in the more detailed curves are separated by 10 3 iterations. 

We note that there is no difficulty to extend the analysis by including the 
complete 0(d 2 ) order in the gradient expansion. The preliminary numerical 
results show no qualitative modification of the behavior presented above. The 
divergence of the coupling constants defined by expanding and -Zfc($) 

increases with the order of $ what suggests that the relevant operator of this 
scaling regime might be nonlocal. 

We close this Section by few remarks about the interior, unstable region. 
The problem in this region is related to the appearance of the saddle points 
at the blocking transformation what correspond to plane waves and make the 
renormalization group step very involved. Apart from the technical complica- 
tions due to the inhomogeneous background field a conceptual problem arises, 
namely the renormalization at the tree level. This happens when there is a 
saddle point with nonvanishing length scaletl what gives an 0(h) contribution 
to the blocking relations. The novel feature of the tree level renormalization is 
that it does not fit into the usual classification scheme of the coupling constants 
what is based on the inspection of the loop integrals. The classical differential 
equations may produce richer and more singular dependence in the coupling 
constants than the polynomes of the perturbation expansion. This actually 
happens in the unstable region and the resulting blocking relation offers no 
guarantee that change in the blocked action will be small, infinitesimal if the 
the cutoff is decreased in a small, infinitesimal amount. By assuming that 
the renormalized trajectory, i.e. the path integral possesses finite derivative 
with respect the cutoff within the unstable region we have an additional con- 
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sistency relations. It simplifies the problem and the important saddle point 
contributions to the blocking can be resummed with the results 

U k ($) = -h 2 <S> 2 + c(k), (32) 

where c(k) is chosen to have continuous potential at the singular curve k 2 = 
k 2 r (Q). Observe that this potential is indeed a "fixed point", i.e. stays invari- 
ant under blocking. For each plane wave saddle point there is a zero mode 
what is related to the translation in the direction of the wave vector. The 
integration over the zero modes finally restores the translation invariance. 

Another, independent support of this form comes from a simple exact 
relation given in terms of the variables in their original, dimensional forncj, 

dV k {$) cW fc ($) _ dW k {$) dv k ($>) , , 

dk d$ dk d$ ' ( ' 

where 

W t (4) = k' + U'^t). (34) 

It was derived by using the gradient expansion only, without making reference 
to the loop expansion. This is a higher order differential equation which has 



spurious solutions but it is at least satisfied by (32). We believe that the 
simplicity of our result, (^), originates from the same balance between the 
energy and entropy as in the case of the Maxwell construction except that the 
role of the domains is played by the zero modes of each saddle point plane 
wave with k < fc cr (0). 



5 Finite Temperature 

The vacuum of the symmetry broken theory is in the stable region so the 
singularities of the coupling constants studied in this work are not obviously 
important for the fluctuations around the vacuum. But the singularity reap- 
pear when the symmetry is broken by controlling external environment vari- 
ables, such as the temperature. Imagine the cooling of a ferromagnet slightly 
above the Curie point or the hot Universe before the spontaneous breakdown 
of a global symmetry group. The order parameter is at its symmetrical value 
<I> = in the high temperature phase and the fluctuations are characterized 
by the coupling constants at $ = what diverge and reflect the nonuniversal 
behavior investigated above. When the system arrives at the symmetry broken 
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phase then the order parameter remains close to the minimum of the effective 
potential what is outside of the spinodial unstable region so long as the time 
evolution is adiabatic. When T « T cr the vacuum is far from the unstable 
region what is detected by the large amplitude fluctuations only. Thus the sad- 
dle points of the inner, unstable region influence the evolution around T « T cr , 
new scaling law can be observed in the vicinity of the phase transition. 

To outline the procedure at finite temperature we note that the blocking 
should be made in three space only in order to preserve the value of the physical 
temperature on the renormalized trajectory. Such an anisotrop blocking in the 
Euclidean space-time where the cutoff k refers to the three space and 
fluctuations in time are eliminated at the one loop level yields the equatior 

k m Um = S logR + fc 2 + ^($)] (35) 



Til-ik 3 f 1 



2(2tt) 3 [T 



;v /A: 2 + + 2 log 



1 



=-t\A 2 +^'(*) 



where uj n = 2%nT. Another simpler strategy is to eliminate the \n\ > jip 
Matsubara modes perturbatively and to apply the blocking for \n\ < uq onrjEil. 
The beta functions derived by either method diverge in the vacuum at $ = 
k = as T — > T cr from above and one expects a singular structure to appear 
what is similar to what was presented above in four dimensions at k — k cr (0). 



6 Summary 

We presented evidences that radically different scaling laws are present in the 
<j) model with spontaneous symmetry breaking at high and low energies. The 
system appears nonuniversal at the phase transition where the influence of 
the nonrenormalizable coupling constants what is suppressed at the ultraviolet 
scaling regime can be compensated for by the singular tree level structure of 
the condensate formation, ft was mentioned that the renormalization group 
method can successfully be applied to the spinodial instability, the mixed 
phase. The saddle points of the blocking procedure which are plane waves 
can be taken into account and the integration over the zero modes restores 
the homogeneity of the mixed phase, what is the reminiscent of the Maxwell 
construction. 

ft is worthwhile noting that as one increases the number of terms retained 
in the potential then the qualitative behavior changes completely when the 
potential was retained up to 0(<I> 22 ). We believe that the divergences occurring 
in the gradient expansion suggests the presence of nonlocal relevant operators 
at the low energy scaling regime. 
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There are several questions left open by these results. What we find the 
most pressing is the classification of the possible nonrenormalizable parameters 
what influences the seemingly nonuniversal dynamics of the phase transition. 
A related issue is the more useful application of the amplification mechanism of 
the divergences generated by the instability as a " renormalization group micro- 
scope" in the coupling constant space to discover the microscopic parameters 
from the long distance observables. 
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